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ABSTRACT 

In this Letter, the properties of the extended radio emission form SDSS DR7 
quasars with 0.4 < z < 0.8 is studied. This low redshift sample is useful since 
any corresponding FIRST radio observations are sensitive enough to detect ex- 
tended flux in even the weakest FR II radio sources. In the sample, 2.7% of the 
sources have detectable extended emission on larger than galactic scales (> 20 - 
30 kpc). The frequency of quasars with FR II level extended radio emission is 
~ 2.3% and > 0.4% of quasars have FR I level extended radio emission. The 
lower limit simply reflects the flux density limit of the survey. The distribution 
of the long term time averaged jet powers of these quasars, Q, has a broad peak 
~ 3 x 10 44 ergs/sec that turns over below below 10 44 ergs/sec and sources above 
10 45 ergs/sec are extremely rare. It is found that the correlation between the 
bolometric (total thermal) luminosity of the accretion flow, L hoh and Q is not 
strong. The correlation of Q with narrow line luminosity is stronger than the 
correlation with broad line luminosity and the continuum luminosity. It is there- 
fore concluded that previous interpretations of correlations of Q with narrow line 
strengths in radio galaxies as a direct correlation of jet power and accretion power 
have been overstated. It is explained why this interpretation mistakenly over- 
looks the sizeable fraction of sources with weak accretion luminosity and powerful 
jets discovered by Ogle et al (2006). 

Subject headings: Black hole physics — galaxies: jets — galaxies: active — accre- 
tion, accretion disks 
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1. Introduction 

Three major questions regarding relativistic jets in radio loud extragalactic radio sources 
that are of fundamental interest are, "how powerful is the typical radio jet?", "how frequently 
do quasars produce long term jets?" and "how is the jet power related to the accretion 
flow thermal luminosity?". The signature of long term powerful jets are the radio lobes of 
Fanaroff- Riley II (FR II) radio sources. Quasars with double radio lobes are conspicuous, 
but quite rare. Thus, radio selected samples are typically used to address the questions above 
since sufficient sample size is required for statistical merit. Most radio l oud sources E do not 



display radio lobes with the resolution and dynamic range of FIRST (IdeVries et al.ll2006l ). 
Thus, in order to find the frequency of the FR II phenomenon in quasars and the complete 
distribution of their kinetic luminosity, one must have the ability to detect lobe emission 
that is below the FR I/FR II divide. The FR I/FR II distinction is both a morphological 
and luminosity classification. In ter ms of long term average jet power, Q, this approximate 



boundary is at Q ~ 5 x 10 43 ergs/sec (IRawlings and Saunderdll99ll ; |Punslyll2008l ). This divide 



corresponds to a 1.4 GHz flux density < 12.5 mJy at z < 0.8, well within the sensitivity 
of the FIRST radio survey, even if the emissivity is distributed over 10 beamwidths (see 
equation (1), below). Thus, the extant databases are now at a point where the questions 
above can be addressed rigorously. 

We intend to study Q with a carefully selected, large sample of SDSS DR7 quasars 
(0.4 < z < 0.8) that is cross-referenced to the FIRST database. Selecting the sources in 
this way reveals prominent narrow lines ([OIII] A5007 and [Oil] A3727) and two broad lines 
H/3 A4861 and Mg II A2798 in the SDSS spectral window for all objects. This allows us 
to compare the various derived properties in the context of a single sample of objects in a 
way that removes the uncertainties associated with sample selection biases if a different set 
of objects and lines are used for the analysis depending on redshift. The cross-correlations 
between the different estimators of the photo-ionizing continuum allow us to segregate out 
the effects of Doppler beaming and line emission created by jet induced shocks, thereby 
reducing systematic uncertainties in the bolometric luminosity, Lf, i, the broadband thermal 
luminosity from IR to X-ray, of the active galactic nucleus (AGN). 

In the next section, the sample is described as well as the extended radio flux data reduc- 
tion and the optical spectra data reduction. The third section summarizes the distribution 
of Q that is found. The fourth section discusses the correlations of Q with the optical/UV 
spectral parameters. 



1 The radio loudness, R, is us ually defined as a 5 GHz flux dens ity 10 times larger than the 4400A flux 
density, R = SsGBx/S ° < 10 jKellermann fc Paulinv-Tothlll969h . 
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Fig. 1. — The histogram of the quasar long term time averaged jet power, Q, computed by 
means of equations (1) and (2). 
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2. Sample Description and Data Extraction Methods 



In order to determine the dependence of line luminosity on the thermal continuum in 
a quasar, we constructed a sample of SDSS DR7 quasars with 0.4 < z < 0.8 and spectra 
with a median S/N > 7, this yielded 10069 AGNs . This base sample has been analyzed 
previously for other purposes in iPunsly and Zhang I (120101 1201 if ). The SDSS DR7 data was 
cross-referenced to the FIRST database. Of these, there were 9270 sources found within 
these FIRST fields that comprise our sample. Most of the radio loud sources did not have 
extended radio flux with the dynamic range and resolution of FIRST. Indicating a sub- 
galactic compact st ructure or very weak , FR I level, extended emission. The optical spectra 
were reduced as in IPunsly and Zhang I ( 120111 ) yielding the following for every quasar: the 
optical/UV continuum luminosity from 5100 A to 3000 A, L cont , the luminosity of the broad 
components of H/3 A4861, Lap, and Mg II A2798, L,Mgii, and the luminosity of the narrow 
lines [OIII] A5007, Loin, and [Oil] A3727, L u- 

The spatial resolution associated with the interferometer beamwidth of the FIRST sur- 
vey is ~ 20 kpc at z = 0.4 and ~ 30 kpc at z = 0.8. The optically thin extended flux density 
on scales of this size a nd larger can be us ed to estimate the long term time averaged kinetic 
energy flux in a jet, Q (IWillott et allll999l ). FIRST observations are at a fairly low frequency, 
1.4 GHz, so they are sensitive to steep spectru m optically thin emission. Resolved emission 
at 1. 4 GHz is very likely to be optically thin (lAntonucci and Ulvestadl Il985t Murphy et al 
19931 ). Thusly motivated, we want to extract this extended flux from the FIRST radio 
database for the SDSS quasars. 

The first step was to cross-reference all FIRST sources within 30" of the SDSS sources. 
Each one of these radio fields was inspected visually. Diffuse emission that is roughly con- 
centric with either the optical quasar position, or an unresolved radio core located within 
3" of the optical position, is considered extended emission of the quasar. This would be 
core - halo type emission. Components within 10" of the optical/radio core are considered 
extended emission. Distant components located 10" to 70" from the quasar/radio core posi- 
tion are considered to be associated with the quasar if there is an extension or jet pointing 
towards the direction of the optical quasar/radio core position, or if there is a second distant 
component located on a line that passes close to this optical quasar/core position (i.e., a 
classical triple). We also considered partially resolved sources in the image plane, quantifying 
the extended emission as the integrated flux minus the peak flux. Other configurations were 
not considered extended flux associated with the quasar (either chance background sources 
or core flux). 

The most accurate estimates of Q should use an isotropic estimat or such as th e radio 
lobe flux. The sophisticated calculation of the jet kinetic luminosity in IWillott et all (119991 ) 
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incorporates deviations from the overly simplified minimum energy estimates into a mul- 
tiplicative factor / that represents the small departures from minimum energy, geometric 
effects, filling factors, protonic contributions a nd low frequency cutof f. The quantity, /, is 
argued to be constrained between 1 and 20. In lBlundell and Rawlingj (120001 ). it was further 
determ i ned th at / is most likely in the range of 10 to 20. T hus, choosing a value of / = 15, 
Punslyj (120051 ) converted the analysis of IWillott et all ( 119991 ) to the formula in equations (1) 
and (2), even though it is just a time average. It is assumed throughout this Letter that 
H o =70 km/s/Mpc, — 0.7 and Q m = 0.3. This formula is an isotropic method that allows 
one to convert 151 MHz flux densities, F 15 i (measured in Jy), into estimates of Q (measured 
in ergs/s): 



Q « 1.1 x 10 45 [(1 + zf +a Z 2 F 151 ] 1 ergs/s , 
Z = 3.31 - (3.65) 

x ([(1 + zf - 0.203(1 + zf + 0.749(1 + zf + 0.444(1 + z) + 0.205] 



-0.125 



(1) 



(2) 



where F 151 is the total optically thin flux density from the lobes (i.e., no contribution 
from a Doppler boosted jet or the radio core). The spectral index is defined by the 
convention, F v ~ v~ a . The appropriate application of this equation requires that one must 
extricate the diffuse lobe emission from the Doppler boosted core and jet. The expression, 
(1), requires 151 MHz flux densities, so we extrapolate the 1.4 GHz extended flux from 
the FIRST observatio n. Without further knowl edge, for lobe emission a ~ 0.8 is the most 
typical value found in iKellermann et al I (119691 ). and this value is used to extrapolate the 
1.4 GHz flux density to 151 MHz. These extrapolations were also cross-referenced to any 
available low frequency data in NED for a consistency check. There is no evidence that a 
significant amount of flux was resolved out by the FIRST beam s ize. For a discu s sion o f 
the errors implicit in equatio n s (1) and (2) , see t he discussions of iFernandes et all (120 10); 
Blundell and Rawlingsl (120001 ); IWillott et all (119991 ) . The independent derivation in iPunsly 
(120051 ) indicates that most estimates in a large sample will be accurate to within a factor of 
2. 



3. Long Term Time Averaged Jet Power 



Figure 1 is a histogram of Q for the 266 quasars with detected extended radio emission. 
Notice that there is a par tition associated with an approximate FR I/FR II dividing line 
at Q = 5 x 10 43 ergs/sec (IRawlings and Saunders! Il99ll ; IPunslyl 120081 ) . The histogram shows 
that the FIRST sensitivity in this redshift range is capable of detecting FR I level extended 
emission. There are 220 FR II sources above the FR I/FR II divide and 46 sources below the 
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divide. This corresponds to 2.3% of the optically selected quasars have FR II level extended 
emission and > 0.4% have FR I level extended emission. The fraction that are FR I sources 
would certainly increase with deeper radio observations. The fraction that are FR II would 
also likely increase with deeper radio imaging (change some FR I levels to FR II levels) as 
faint diffuse flux is likely swamped by the FIRST rms noise. Th e frequency of 2.3% f or the 
FR II quasars is consistent with the value of 1.7% deduced by IdeVries et al.l (120061 ) since 
their sample has a mean of z ~ 1.3, so many faint, diffuse FR II sources (that are prominent 
in Figure 1) would be difficult or impossible to detect with FIRST sensitivity at these high 
redshifts, hence their lower deduced rate of occurrence. The significant fraction of FR I 
extended fl ux densities is consistent with the discovery of quasars with F R I extended radio 



emission m 



Gower and Hutchingsl ( 119841 ); iBlundell and Rawlings I ( 120011 ) . 



A major advantage of this study is that the peak of the quasars extended luminosity 
distribution is resolved in Figure 1 due to the FIRST sensitivity in this redshift range. It 
shows that the most frequent luminosity (for sources with extended emission) is definitely at 
the FR II level. The peak level of the histogram is > 10 44 ergs/sec with marginal evidence of a 
second peak at ~ 5x 10 44 ergs/sec. Only 33 sources have Q > 10 45 ergs/sec, ~ 0.3%. Although 
powerful FR II quasar s (Q > 10 45 ergs/sec ) are prominent in low frequency radio surveys 
such as the 7C survey, IWillott et ajJl999h . these sources are extremely rare, perhaps the 
rarest subcla ss of quasar - an o rder of magnitude rarer than low ionization broad absorption 
line quasars (jZhang et alll2010l ). Whatever physical phenomenon creates these large values 
of Q must be very difficult to achieve and maintain for long periods of time in a quasar 
environment. 



4. Correlations Between Accretion Luminosity and Jet Power 



The blue continuum luminosity is the most basic signature of the thermal emission 
from the quasar, so it is the most commonly used quantity for estim ating L^ni.. Perhap s 
the most popular bolometric correction is the simple one proposed by iKaspi et all ( 120001 ). 
Lboi ~ 9ALa(5100 A). Clearly, using a portion of the optical/U V continuum is mor e accu rate 
than a single point and we have that at our dispo sal, L rnT1t . In iPunslv and Zhang I (120111 ) , in 
another study of these 10,069 SDSS sources, the Kaspi et al (12000 ) relation was converted 
to 

Lboi ~ 15L cont . (3) 

L cont is the most direct way to estimate Ly. One can also estimate Lboi from the spectral 
luminosity at 3000 A for the sources in this sample. The average spectral index of the 
optical/UV continuum for the radio quiet sources in the sample of 10,069 quasars was found 
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in iPunsly and Zhang I (j201lf ) to be a 

estimator can also be realized as 



0.55, where L h 



nu 1/ 



^bol 



7.1AL A (3000A) . 



Thus, the Kaspi et all fl200Q[ ) 

(4) 



It is interesting to see i f there is a strong relationship between L^g and Q as has been 
inferred for radio galaxies in iRawlings and Saunders! (1l99l[ ) ; IWillott et all ( 119991 ) . An uncer- 
tainty in equation (3), in the the context of radio sources, is the high frequency tail of the 
synchrotron emitting jet and its potential to dilute the pure L cont from the accretion flow. 
Our carefully selected sample ameliorates these concerns. We have other reliable surrogates 
for Lb i. The synchrotron tail dies of rapidly with frequency in the opti cal/UV and the big 
blue bump from the accretion disk becomes more prominent in the UV (IMalkan and Moore 
19861 ). The spectral index of L u of the optical/UV synchrotron tail in quasars is typically, 
a ~ 2, compared to the a ccretion spectral luminosity with a ^ 0.55 f or the 10,069 sources in 
the SDSS quasar sample ( jSmith et allll988l ; IPunsly and Zhang 11201 ll ). Thus, the ratio of the 
synchrotron tail luminosity (XL\) to the accretion luminosity is typically > 2 times larger 
at 5100 A than at 3000 A. Consequently, if synchrotron dilution is significant, there should 
be a noticeable difference in the correlati ons of quasar parameters w ith L cont and those with 
L(3000A). Furthermore, it was shown in IPunsly and Zhang I ( 1201 ll ) that the broad compo- 
nents of Mg II and in particular H/3 are excellent surrogates for the underlying accretion disk 
luminosity. These spectral parameters are available for all sources in our sample. Figure 2 
shows that L H/3 is an excellent surrogate for L cont for our radio sources as well. The best 
fit relationship displayed at the top of the plot is nearly linear. The tight correlation is 
evidence that Doppler beaming of the high frequency tail of the synchrotron jet is not a 
major contaminant to the continuum luminosity measured by SDSS in a statistical sense. 

We explore the correlations between the spectral features and Q in the correlation 
matrix, Table 1. The statistical metric of choice is the Spearman rank correlation coefficient 
in Table 1. As expected, L cont and L (3000 A) are tightly correlated and their correlations 
with Lh/3 are also tight (see Figure 2). The most intriguing aspect of Table 1 and Figure 
3 is that the weakest correlation with Q is from the continuum itself and then Lnp- The 
correlations of Q with the narrow line strengths are more significant than with the broad 
lines or the continuum. A lack of a s trong correlation betwe en Q and L cont is consistent with 
the right hand panel of Figure 4 of iFernandes et all (120101 ) based on a sample of optically 
selected radio loud quasars. 

Many previous stud ies have shown a strong correlation between narrow line strengths 



and Q for radio galaxies ( 


Rawlines and Saunder 


s 1991 


2010. 


2011; 


Koziel-Wierzbowska and Stasinska 


2011). 



of the continuum in Table 1 (the broad lines and a direct view of the continuum itself) 
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Table 1: Spearman Rank Correlation Matrix for the Quasars with Extended Emission 

-^cont £(3000 A) Lh/3 L/Mgll Loin Lou Q 

1 0.985 0.869 0.732 0.525 0.379 0.314 

L(3000A) 0.985 1 0.871 0.775 0.522 0.367 0.313 

L Hp 0.869 0.871 1 0.799 0.636 0.468 0.417 

L Mg ii 0.732 0.775 0.799 1 0.565 0.397 0.421 

Loin 0.525 0.522 0.636 0.565 1 0.844 0.507 

Lou 0.379 0.367 0.468 0.397 0.844 1 0.467 

Q 0.314 0.313 0.417 0.421 0.507 0.467 1 



Hp LINE STRENGTH VERSUS THE CONTINNUM LUMINOSITY 

LOG[Hp Luminosity (ergs/sec)] = 1.0134{LOG[Continuum Luminosity (ergs/sec)]} - 2.152 
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Fig. 2. — The scatter plot of log Lap versus logL cont for the 266 quasars with extended 
luminosity. The dashed lines are the 2.5 a level scatter from the linear fit that were chosen 
to highlight the small number of sources below the trend-line that might be significantly 
affected by Doppler boosting. 
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eit her did not exist o r were hidden by the dusty torus ( lAntonuccil Il993l ) . It was shown 



m iPunsly and Zhang I (120111 ) that a narrow line luminosity is not nearly as accurate of a 
measure of L bo i, as L cont or a broad line luminosity. It seems that t he interpreta t ion o f 
the correlation of L n with Q for the relatively powerful 7C sources in IWillott et all (119991 ) 
in terms of a strong correlation between and Q might need some revision. In Table 
1, the correlations of L OII and L OIII are stronger with Q than with L cont or L (3000 A). 
This suggests that jet propagation is a comparable or more significant source of excitation 
energy than photo-ionization for the O II emitting gas in FR II radio sources. It has been 
shown that Loin is superior to Lou as a measure o f Lboi , due to at least two effects: 
the dependence on the io nization parameter. ISimpsonl (119981 ). and contamination by star 
formation in host galaxies (IHeckman et alll2004l ). Thus, correlations with Loin have a more 
direct physical interpretation. For O III, the source of excitation for the narrow line gas 
might be roughly an equal mix of photo-ionization energy and jet induced energy based on 
the correlation strengths in Table 1. Thus, the narrow line strength and Q are dependent 
statistical quantities. So, if a narrow line strength is used as a surrogate for Lboi in a 
correlation analysis or scatter plot with Q then the probability of large outlier sources (i.e., 
large Q and small narrow line strength) is greatly diminished compared to a scatter plot in 
the Q - Lboi plane (i.e., large Q and small Lboi) because Q is itself a strong source of narrow 
line emitting gas. This interpretation seems to rectify a discrepancy seen here and elsewhere 
in the lit erature. Consider the bimodal scatter in the IR luminosity/ radio luminosity plane 
noted by I Ogle et all (120061 ) in their Figure 1. They concluded that significant Q can exist in 
sources where there is essentially very little accretion thermal luminosity as well as in sources 
with large (quasar-like) Lboi. Conv e rsely, the cor r elatio n between narrow line strength and 
Q in iRawlings and Saunders! (Il99ll ) ; IWillott et all ( 119991 ) led to the conclusion that accretion 
power was very closely linked to jet power. Based on Table 1, this difference in interpretation 
arises because a strong radio source with a weak ionizing continuum is still likely to be a 
major source of narrow line emission, thereby keeping these sources correlated with narrow 
line strength - even though they are not correlated with L^oi- An important point to bear in 
mind is that only a small range of optical luminosity is probed here, with most sources being 
within one order of magnitude. Given the very large scatter in some of these correlations, it 
is clear that over such a small range of luminosity, correlation statistics can give low values. 
Observations over many orders of magnitude with the same scatter would giv e much higher 



correl ation coefficients, e.g. the correlations of IRawlings and Saunders! ( 1 199 if) ; IWillo 



(119991 ) . However, the scatter in these relationships is still significant, 0.54 dex (IWillott 



fetal 



2000). 



- 10 - 



JET POWER AS A FUNCTION OF CONTINUUM LUMINOSITY 
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Fig. 3. — On the left of the first row is the scatter plot of logQ versus logL cont and on the 
right is the same for log AL A (3000 A) for the 266 quasars with detected extended luminosity. 
Since the two scatter plots look very similar, the effects of Doppler beaming must be negligible 
within the sample as discussed in the text. Scatter plots of the logarithm of the line strengths 
as a function of the logarithm of Q are in the next two rows. The narrow lines show a modestly 
elevated level of correlation relative to the broad lines in spite of some outliers with low line 
luminosity. The largest scatter is with the continuum luminosity. 
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5. Conclusion 



This Letter used a large SDSS DR7 sample and FIRST radio observations to address 
three fundamental questions in the study of quasar radio jets. We repeat them here with 
answers indicated by the sample considered: 



How powerful is the typical radio jet? From Figure 1, typically, the long term 
time averaged jet power, Q, is a factor of a few above the FR I/FR II divide, Q ~ 
1 — 5 x 10 44 ergs/sec. 

How frequently do quasars produce long term jets? For 0.4 < z < 0.8, 2.3% 
of the optically selected quasars have FR II level extended emission and > 0.4% have 
FR I level extended emission indicative of long term jet production. 

How is the jet power related to the accretion flow thermal luminosity? In 

the last section (see Figure 3), it was established that the long term time average jet 
power, Q, is not strongly correlated with the accretion flow thermal luminosity. The 
extended lobe emission used to compute Q was emitted from the quasar central engine 
10 5 - 10 7 years earlier than the observed optical/UV continuum used to evaluate 



( iBlundell and Rawlingsl 120001 1 . Thus, the two properties are not contemporaneously 
determined and there is no apriori reason to expect them to correlate strongly, except 
possibly by the overall scaling of the central black hole mass responsible for both. 
No information was found in the sample that can be used to reliably determine a 
connection between the instantaneous jet power and Lboi- 



Another major finding of this work, that is closely related to point 3 above, is that 
powerful FR II quasars (Q > 10 45 ergs/sec) are extremely rare, ~ 0.3% of all quasars. Thus, 
any theoretical model of the central engine must explain why it is so difficult for a quasar 
environment to maintain a configuration conducive to strong jet formation over ~ 10 5 - 10 7 
years. 



We would like to thank the referee for many cogent and useful comments. 
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